Universal Correlations and Dynamic Disorder in a Nonlinear Periodic ID System 
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When a periodic ID system described by a tight-binding model is uniformly initialized with equal 
amplitudes at all sites, yet with completely random phases, it evolves into a thermal distribution 
with no spatial correlations. However, when the system is nonlinear, correlations are spontaneously 
formed. We find that for strong nonlinearities, the intensity histograms approach a narrow Gaussian 
distributed around their mean and phase correlations are formed between neighboring sites. Sites 
tend to be out-of-phase for a positive nonlinearity and in-phase for a negative one. The field 
correlations take a universal shape independent of parameters. This nonlinear evolution produces 
an effectively dynamically disordered potential which exhibits interesting diffusive behavior. 



The tight binding approximation is one of the simplest 
models that predict band structure and ballistic motion 
in periodic systems Using a tight-binding model with 
disorder, P.W. Anderson was able to predict and study 
the effect of localization in disordered lattices 0|. A non- 
linear version of the tight-binding model, better known as 
the Discrete Nonlinear Schrodinger Equation (DNLSE) 
has been used to study nonlinear evolution in periodic 
systems, initially in the context of periodic molecular and 
mechanical systems and extensively in recent years 
to describe nonlinear propagation in optical waveguide 
lattices 13, 0|, as well as matter- waves in light-induced 
lattices(6|]. In particular, the DNLSE explains the forma- 
tion of nonlinear intrinsic localized modes, also known 
as discrete solitons 0| or discrete breathers Q. Disor- 
der and localization have been studied experimentally 
recently both in matter-waves [j| [13] and in photonic 
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f particular interest and some 
debate is the interplay between localization and nonlin- 
earity dSEii. 

Here we report on a new phenomenon that results from 
the interplay of nonlinearity and disorder. This effect oc- 
curs in periodic systems, where nonlinearity induces dis- 
order in an otherwise perfectly ordered lattice. What we 
find is that when the system is initialized with random- 
phase fields, it evolves into particular distributions with 
well defined stationary statistical properties. Most in- 
terestingly, the field correlation function and the distri- 
bution of phases assume universal forms independent of 
the exact value of the nonlinear parameter. The result- 
ing distribution induces a dynamic structure with several 
intriguing properties. 

In Fig. 1 we show the results of an optical experiment 
in a waveguide array that motivated this study. Light 
with uniform intensity yet random phases was injected 
into a large number of waveguides in a periodic waveg- 
uide lattice. The intensity at the output end was mea- 
sured, and the histograms of intensity values obtained 
from many repeats of the experiment with different ran- 
dom phase realizations are shown for both low inten- 
sity and high intensity light. For experimental details, 




Figure 1: Experimental measurement of output light inten- 
sities from a waveguide array, when input fields with equal 
amplitudes, yet random phases, are coupled to several ad- 
jacent waveguides. Measured values of many random-phase 
realizations are shown as histograms for linear propagation 
(blue circles) and high-intensity, nonlinear propagation (red 
circles). The inset shows two sets of measurements for an 
input with the same phase realization after linear (blue) and 
nonlinear (red) propagation. The output intensity distribu- 
tion is exponential-like in linear propagation, and becomes 
narrow, Gaussian-like in the nonlinear case. 



see Methods below. While linear propagation produced 
exponential-like distribution of intensities, as expected 
for summing of many random-phase inputs, the nonlinear 
propagation yielded a much narrower distribution around 
the average intensity. Increasing the optical power led to 
the narrowing of the output distribution. Note that for 
both measurements, the counts at low intensity values 
are underestimated because of scattered light and instru- 
mentation noise. To investigate this behavior we have 
modeled the problem using the DNLSE. While we use 
here the notations of optics, our results will hold in gen- 
eral for all other systems described by the DNLSE. The 
evolution of light in a periodic array of weakly coupled 
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waveguides is described by: 

i^H. = C(a n _i + a„+i) +7|a„| 2 a„, (1) 

where a„ is the amplitude of the mode in the n th waveg- 
uide, C is the coupling coefficient to the nearest neigh- 
bors and 7 is the nonlinear coefficient, positive (negative) 
for focusing (defocusing) nonlinearity. We shall consider 
the situation where light is injected into the array with 
uniform amplitudes \a n \ = ao, yet with completely mi- 
correlated, random phases. It is convenient to use the 
normalized equation, 

i— if- = {u n -i +u n+1 ) +T\u n \ u n , (2) 

where ( = zC, u n — a n /ao and T — jao 2 /C. With this 
normalization, the input intensities are all uniform with 
7 n (0) = u n u n * = 1. 

Consider first the linear problem, i.e. r = 0. As might 
be expected, after a certain distance, mixing of the dif- 
ferent input fields leads to an output pattern with fluctu- 
ating intensities. Fig. 2(a)-(c) shows results of numerical 
simulations of Eq. (2) for various properties of the fields 
after propagating a distance of C = 10 in an array with 
TV = 256 waveguides. Periodic boundary conditions are 
used to avoid edge effects. The results shown are aver- 
aged over 500 realizations with different random initial 
fields. Fig. 2(a) shows the intensity histogram; it fol- 
lows an exponential law, P(I) = cxp(— I), as expected 
for random fields. Fig. 1(b) shows the field correlations 
nU n j r ] < . ± u n u n +k) /2N demonstrating, as could 
be expected, that the fields at different sites do not cor- 
relate. Finally, Fig 1(c) shows the histogram of phase 
differences between neighbors, 9 n = <p n — (f> n +i, with 
4> n = argu„ the phase of the field u n . These phases 
are uniformly distributed. 

We now repeat the simulations with nonlinearity, and 
we will be interested mostly in the limit of strong nonlin- 
earity; results for T = ±20 and for T = ±200 are given 
in Figs. 2(d)-(f). Two changes from the linear case are 
obvious. First, the intensity histograms, shown in Fig. 
2(d), now converge around the average intensity value 
of 1, with a width that shrinks with the nonlinear pa- 
rameter. The distribution seems to fit well a gaussian 
distribution with P(I) = exp[— (I - l) 2 /2a 2 ), and it is 
independent of the sign of the nonlinearity. 

The second major effect of the nonlinearity is the in- 
duced spatial field correlations. Most interestingly, the 
correlation function (Fig. 2(e)) takes a shape that is in- 
dependent of the nonlinearity value, and is sensitive only 
to its sign. It is described well by exponential decays, 
Ck = (— l) fc cxp(— ak) for positive (focusing) nonlinear- 
ity and Ck = exp(— ak) for the negative case. Note that 
the correlation is only visible in the fields - the intensi- 
ties remain uncorrelated; Intensity correlations show a di- 
minished peak at k=0 and a uniform background. These 




Figure 2: Simulation results of the DNLSE with uniform in- 
tensities and random phases, (a) Intensity histogram, (b) field 
correlation and (c) phase-difference histogram for a linear sys- 
tem (r = 0), exhibiting exponential intensity distribution and 
uncorrelated fields and phases, (d) Intensity histograms for 
r = ±20 (blue circles) and F = ±200 (red circles). The distri- 
butions are narrower at higher nonlinearities, yet independent 
of the sign of V. (e) The universal field correlation function 
is identical for both nonlinear values, but depends on their 
sign: blue for positive (focusing, T > 0), red for negative 
nonlinearities. (f) Phase difference histograms also approach 
a universal distribution, concentrating around n for positive 
nonlinearity and for the negative case. The theoretical lines 
in (d)-(f) are predictions of Eqs. (7) and (8). 



field patterns are consistent with the known properties of 
modulation instability in such systems: Staggered (un- 
staggered) fields are stable in positive (negative) nonlin- 
earity arrays [151 ]. 

Since the intensities become more uniform at high non- 
linearities, the field correlation functions are dictated by 
the variations of phases between neighboring sites. In 
Fig. 1(f) we show the histograms of these phase differ- 
ences for positive and negative nonlinearities. While in 
the former neighboring waveguides are most likely to be 
out of phase, as the distribution peaks at n, in the latter 
two neighboring waveguides tend to be in-phase. This 
distribution of phases also attains a constant profile at 
high nonlinearity values. 

The field correlation and the intensity distribution are 
closely related. This can be deduced from the conserved 
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quantities, the Hamiltonian 

h = ir^/ 2 + J]K< +1 + < u „ +1 ), (3) 

and the total photon number, 

A = J2 J n- (4) 

Prom these it is easy to show that 

jo 2 (Q+C 1 (Q = H (5) 

is also constant. Here a — (El%/N — l) 1 / 2 is the standard 
deviation of the intensities. Since our initial condition are 
of uniform intensities, er 2 (0) = 0, and random phases, 
Ci(0) « 0, then H = 0, hence 



^ 2 (0 = -Ci(C). 



(6) 



Equation (6) predicts that the signs of C\ and T are dif- 
ferent, as indeed is observed in Fig 2. For weak non- 
linearity (r <C 1), when the distribution deviate only 
slightly from exponential, a small correlation is formed 
with C\ = — r/4. However, as the nonlinearity increases, 
Eq. J6J predicts that the intensity distribution has to 
narrow down, since necessarily \C\\ < 1. 

To relate the phase and intensity fluctuations, we have 
to study the statistical properties of the DNLSE. Such 
investigations were carried out by Rasmussen et al 
and extended later by Rumpf [17|. They were mostly in- 
terested in the conditions for the generation of localized 
structures. With our initial conditions, it can be shown 
that localized structures are not formed, but we can use 
the same formalism to derive the probability distribu- 
tions for phase and intensities. 

In essence, the state of maximal entropy 
S\p(Ii, ...In, #1, ■■■On)] = — J pln.pY[dIid9i can be de- 
rived by the variational problem S(S — aA— r\H — A J p) = 
0, where a,rj and A are the appropriate Lagrange mul- 
tipliers The problem is greatly simplified by the 
approximation E(m„u* +1 + c.c.) « 2^cos(6*„), which is 
consistent with the observation that the intensities are 
uncorrelated and at high nonlinearities they are close 
to their average value of 1. With this approximation, 
the intensities and phases are separable, and their 
distributions are derived to be: 



Al exp[-^(/ 



Pl(I) - A, <<xp —!/ - 1) 
pe(0) = A2 exp[-2?7cos(0)] 



(7) 
(8) 



with Ai, A2 appropriate normalization constants and r\ « 
±0.533 is the solution of 477 J cos(9)pg(9)d9+l = 0, where 
the sign is selected to match the sign of T. The phase 
distribution is then maximized at 9 = it (9 = 0) for posi- 
tive (negative) T, respectively. Note that these universal 




Figure 3: Maps of nonlinear potential induced by the fluc- 
tuating fields for (a) T = 20 and (b) F = 200 and the cor- 
responding correlation maps Y(k,5() (c,d). The maps show 
a section of 40 waveguides (vertically) and a propagation of 
C = 1- Note the faster dynamics for the higher nonlinearity, 
evident also in the correlation map. For easy viewing, the 
potentials are normalized to their peaks, although the values 
in (b) are about 3 times higher than in (a) 



phase distribution functions, which are shown as lines in 
Fig 2(f), are independent of the value of the nonlinear- 
ity. They lead to the universal correlation function that 
decays with Ck/Ck+i — — 4ry as shown in Fig 2(e). 

Finally, we would like to discuss the potential (or in- 
dex) pattern which is induced by the fluctuating fields 
through the nonlinear effect. The fields induce an ef- 
fective dynamic disorder into the lattice; While we have 
shown above that the variations in intensities tend to 
diminish at high nonlinearities [a — (^r) -1 / 2 ], the vari- 
ations in induced nonlinear potential, i.e. Ta actually 
increase. The array becomes effectively more disordered: 
A stationary structure with this level of disorder would 
have been characterized by a localization length that 
is narrower than the lattice spacing. However, this in- 
crease in disorder is also accompanied by faster dynam- 
ics: The induced pattern changes faster at higher non- 
linearities. The nonlinear potential map TI n (Q is shown 
in Fig. 3, together with the intensity correlation maps 
Y(k,6() = E n / d(I n (()I n - k {( - 5(). Note that while 
at a given ( the intensities at adjacent waveguide do not 
correlate, Y(k,0) = Sk,o, they are correlated at other 
values of £• It is this dynamic potential structure that 
determines the field correlations and other peculiarities 
of this system. 

One way to characterize this structure is to investigate 
the transport properties of a a secondary probe beam 
that propagates in the nonlinearly induced disordered 
potential. We have simulated the propagation of such a 
weak probe field in an array with 1024 sites. The probe 
is launched into a single central site at £ = 0, and its 
width is monitored along propagation. Fig. 4(a) shows 
this width, averaged over several realizations, as a func- 
tion of C 1 ^ 2 , f° r two values of nonlinearity, r = 20 and 
r = 200, and Fig 4(b) shows the probe intensity distri- 
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Figure 4: The diffusion of a weak probe in the nonlinear po- 
tential induced by the fluctuating fields with F = 20 and 
F = 200. (a) the averaged probe width as a function of C, 1 ^ 2 . 
(b) The averaged probe profile at C = 2000 



bution at £ = 2000. It seems that the broadening is gov- 
erned mostly by diffusive dynamics, that is, a gaussian- 
like distribution that broadens diffusively as C 1 / 2 . In- 
deed, dynamic disorder that is spatially uncorrected is 
known to lead to diffusive broadening [lil [l9| . What we 
find interesting is that at higher nonlinearities the diffu- 
sion coefficient is actually larger, in spite of the stronger 
disorder. This is most likely the result of the faster dy- 
namics, but it could be that the nontrivial correlation 
maps also plays a role. These points are currently under 
investigation. 

In conclusion, we have shown that when systems de- 
scribed by the DNLSE are initialized with equal ampli- 
tudes yet random phases, universal field and phase cor- 
relation are formed in the high nonlinearity limit. In 
contrast with the thermal distribution of intensities ob- 
tained in linear propagation, the intensity variations are 
diminished, and universal phase correlations are formed. 
The underlying potential field leading to this behavior is 
of a dynamic disordered system, which is characterized 
by diffusive propagation. 

Methods: The periodic waveguide lattices were pre- 
pared on an AlGaAs substrate using e-beam lithography, 
followed by reactive ion etching. The width of each of the 
waveguides as well as their spacing were 4 microns. The 
tunneling parameter between sites C was measured to be 
280m -1 . The light source was a pulsed optical paramet- 
ric amplifier, producing 1.2 ps pulses at a wavelength of 



1530nm with a 1 KHz repetition rate. The spatial phase 
of the laser light beam was randomized using a computer- 
controlled spatial light modulator. The beam was shaped 
using a cylindrical optics and then injected into the in- 
put facet of the array, covering 15 lattice sites. The light 
intensity profile was measured at the output facet of the 
8mm long sample using an infrared camera and the in- 
tensity statistics was calculated by collecting many such 
images with random phase realizations. 
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